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Abstract 

A heat conduction problem is studied using extended hydrodynamic equations obtained from Enskog's 
equation for a simple case of two planar systems in contact through a porous wall. One of the systems 
is in equilibrium and the other one in a steady conductive state. The example is used to put to test the 
predictions which has been made with a new thermodynamic formalism. 

PACS: 51.10.+y 05.20.Jj 44.10. +i 05.70.Ln 
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1 Introduction 

Today Enskog's original kinetic theory is known as the standard Enskog theory (SET) 012 El because after 
the pioneer work of van Beijeren and Ernst ^ there are several new versions of Enskog's theory coUectively 
called revised Enskog's theory (RET) 0. Among the latter there are versions that have been extended to 
describe condensed matter 6 . To Navier-Stokes level both SET and RET lead to the same results [7] , 
whether or not an external force is present. 

In the present article we make use of extended hydrodynamic equations for the bidimensional case 
They are more complete than a linear approximation but still they are the result of an approximation scheme 
that we explain elsewhere. Using a strategy as in Ref. 'S' and approximations defined in Rcf. (9|| we obtain 
in Sec. [3 the same hydrodynamic equations for SET and RET. 

In this article we apply our extended hydrodynamics to a one dimensional steady heat conductive state. 
There is much work on this as, for example, the experimental results in or the theoretical ones in jlll 
1121 118j . Recently Kim and Hayakawa studied this problem for hard core and Maxwellian particles 
using Boltzmann's equation combined with Chapman-Enskog's method. They tried a test and criticized the 
analysis of the nonequilibrium steady state thermodynamics (SST) proposed by Sasa and Tasaki jl5| . In 
the last reference the authors state that if there is gas in a one dimensional heat conductive configuration 
in contact, through a porous wall, with an equilibrium gas state, then a pressure difference must appear in 
the direction of the heat flow. We analyze this double system making use of the extended hydrodynamic 
equations derived from Enskog's equation using Grad's moment expansion method (16j. Our conclusions 
differ from those in |15| . 

The organization of the present article is as follows. In Sec. |21the configuration of these systems is drawn 
schematically, in Sec. |31the condition for the two systems to be in contact via the central porous plate is 
introduced: the upper and lower plates are normal plates; the central plate has many small pores through 
which the gas can pass. In Sec. 0]we give the basic equations used in this paper. Comments are in Sec. El 
Finally, our discussion and conclusions are written in Sec. El and (3 respectively. 

*| published in Physica A, 354. 77-87, 2005 
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2 Definition of the system 



Sasa and Tasaki |15| proposed an interesting system consisting of a nonequilibrium steady state subsystem 
in contact with a subsystem in equihbrium as explained below. This system has three plates as shown in 
Fig. n and there is gas between them. The upper and lower plates (plates 1 and 3) are normal plates. The 
central plate (plate 2) has pores through which gas can pass. 

Following Sasa and Tasaki, we consider the system consisting of three infinite parallel plates 1, 2 and 3 
separated by a distance L. The Y axis is defined perpendicular to them while an X axis is placed on plate 
2. The pores in plate 2 are distributed homogeneously. Plates 1 and 2 have fixed temperature Ti while plate 
3 has a different (fixed) temperature T2. 

After a sufficiently long time, by effusion, some of the gas passes through plate 2 and the gas between 
plates 1 and 2 reaches an equilibrium state. The system between the plates 2 and 3 reaches a nonequilibrium 
steady state with translation symmetry along the X axis. 

We assume that the typical distances between pores is very small and that the diameter of the pores is 
also very small, so that the ratio between such lengths and the mean free path is much smaller than unity. 
Having no external force and no hydrodynamic velocity there is no heat flux parallel to the plates. The 
system is in a static configuration. 
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Figure 1: The plates of the system as described in the text. 



3 The contact condition 

In general, there is a difference between the temperatures of the plates and the temperatures of the gas in 
contact with them, this is a well-known effect called thermal slip. However, for simplicity sake, we assume 
that the temperature of plate 2 and the gas in contact on both sides of it are equal, namely we are neglecting 
the Knudsen layer. 

The velocity and the peculiar velocity of the gas will be denoted by c and C, respectively. The condition 
that there is no mass flux through plate 2 is 
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where /oquii. = fieq 2tiTi ^'^P ~W\ equilibrium distribution function associated to the gas between 

plates 1 and 2 and Ueq is the uniform number density in this same region while fy=o is the nonequilibrium 
distribution function between plates 2 and 3 evaluated at ?/ = 0. Next, it is necessary to see how the two 
distributions satisfy condition 
Imposing condition ^ yields 
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In addition, the total mass conservation law for the system is 

n{y)dy ^ 2. (3) 



Above we are using dimensionless fields and dimensionless variables in general. The fields n (number 
density), Pij (pressure tensor), Q (net heat flux vector), and T (temperature) generally depend on the 
coordinate y, where Pij and q*^ are the symmetric and traceless part of the pressure tensor and the kinetic 
part of the heat flux vector respectively. These hydrodynamic fields are defined according to the following 
sum rules: 
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fdc = n{y), (4) 
C{y)fdc = 0, (5) 



lc{yffdc = n{y)T{y), (6) 



Ci{y)Cjiy)fdc = n{y)T{y)dij+Pij{y), (7) 



^C{y)C'{y)fdc = q^^iy). (8) 
We also use the following dimensionless numbers 

Ki = 1. ^ Knudsen number, S = ^ = Kapo 
n L L 

where a is the particle's diameter, £ their mean free at equilibrium and po is the mean area density. 



4 Balance equations 



The basic concrete equations solved here are the following. 

• In the case of the linearized Boltzmann-Grad method (LBG): Pyy{y) = Pyy = constant, Pxy{y) 
Pxy = constant, Qy{y) = Qy — constant, 

n{y)T{y) = P,, = Pyy, 
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(9) 



• In the case of the Enskog-Grad method (EG): Pyy{y) = Pyy = constant, Qyiu) = Qy = constant, 
Pxy{y) =Pxy{y) = 0, 
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The substitution of (5 = and x = 1 in the EG equations yields the equations corresponding to the 
nonlinearized Boltzmann-Grad method (NLBG). 



5 The pressure difference between the equihbrium and nonequi- 
Ubrium sides 

All the results we describe in what follows were obtained using perturbation methods choosing T2 > Ti and 
using e — (T2 — Ti)/Ti as the perturbation parameter. We solve the system of equations and their boundary 
(contact) conditions up to e^. We choose S = 0.001. In such case Kn is in inverse proportion to pQ. We 
choose Henderson's expression ^7] as the concrete expression for x- 
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(14) 



We calculate the pressure in both sides of plate 2. Using Pyy for the nonquilibrium steady state side and 
the pressure Peq which is estimated by the state equation for the equilibrium side 
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(15) 



the pressure difference AP is defined by 



AP — P — P 

LJi.r — ryy j-eq ■ 



(16) 



Note that n(0) = rigq and AP = to first order in e. Hence, we rewrite our results in the following way 



n(0) = rie 



1 + A. 



AP = Aap^ 

where An and Aap are constants. Furthermore, it is possible to rewrite Pyy'. 
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Tables n and 13 give the values of these constants for e — 0.05 and 0.1, respectively. The value and sign 
of AP depend on e and Kn. Table 13 gives the value of Aap obtained by first order EG, that is, up to 5 for 
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An X lU 


/\AP X lu 


\vy 1 \ 

Ap 1 An 


0.005 


-0.360727 


1.43059 


-3.45019 


0.01 


0.6175405 


-0.37938 


-1.17319 


0.02 


0.5891347 


-0.51577 


-1.66426 


0.05 


0.4851223 


-0.46437 


-1.87508 


0.1 


0.4309258 


-0.44057 


-1.93856 


0.2 


0.4158376 


-0.41159 


-1.96948 


values of A„, Aap and X-jy /\n for e 
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Aap X 10-2 


\vy 1 \ 


0.005 


-0.405391 


1.49646 


-3.19368 


0.01 


0.614064 


-0.37163 


-1.16022 


0.02 


0.590175 


-0.51564 


-1.66031 


0.05 


0.486046 


-0.40465 


-1.87365 


0.1 


0.431343 


-0.43093 


-1.93785 


0.2 


0.416123 


-0.41182 


-1.96124 



Table 2: The values of A„, Aap and X^^ /\n for e — 0.1 in the case of EG 



e — 0.05 and 0.1, respectively. In this case, the pressure difference also exists and its value and sign depend 
on e and Kn, too. 

In the case of LBG, since p^x — then A„ = Aap = A^^ 0. There is no pressure difference in this case. 

On the other hand, for the case of NLBG, the substitution of (5 = in Eqs. 101, (EHl, EJ; OH; 03; 
(Cni, (CU and (Uni leads to A„ = 1/256, Aap = -1/256 < and \yy = -1/128 where for A„ and Aap 

it is correct to consider only up to second order in Qy. Hence the osmotic pressure difference does exist and 
its value is constant and negative. 

Furthermore, for the case of EG, A^^/A„ = —2 — S^n^q ^ —2. 

We analyze the pressure difference AP from another point of view. It is sufficient to calculate AP up to 
e^. It is given by 
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It is seen that the sign of AP changes from negative to positive approximately at po = 0.2, whereas it is 
always negative in the NLBG and to first order in the EG's case. 

Besides the system far from equilibrium, we are also interested in a region extremely close to the equi- 
librium condition. Therefore, we analyze the case without the strong nonlinear term, namely, we eliminate 
the terms involving qy{y)^ and Pxx{y)qy[y) in Eqs. ifT^ and (fT^ . In this case, up to (5, 
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Assuming that the derivative dn{y)/dy of the density at plate 2 has the same sign as Qy (this is normaly 
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Kn 


Aap X 10-2 e = 0.05 


Aap X 10-2 e = 0.1 


0.005 


1.208777 


1.269384 


0.01 


-0.417305 


-0.410453 


0.1 


-0.431207 


-0.431631 



Table 3: The values of Aap for e = 0.05 and 0.1 in the case of the first order EG 



correct), AP is always negative. Evaluating up to 5^ and yields 
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As < po < 1, is positive but -f- ^ -2 

Furthermore, we calculate A„, Aap and up to 5'^ and e^. Tables 0] and give the values of these 
constants for e = 0.05 and 0.1, respectively. The value and sign of AP depends on e and Kn, too. 



Kn 


A„ X 10^2 


Aap X 10^2 


A™/A„ 


0.005 


-0.552682 


1.59241 


-2.72540 


0.01 


0.3430882 


-0.16459 


-1.00243 


0.02 


0.2577696 


-0.22015 


-1.63486 


0.05 


0.1182556 


-0.11280 


-1.87030 


0.1 


0.0612808 


-0.05996 


-1.93731 


0.2 


0.0311473 


-0.03082 


-1.96915 



Table 4: The values of A„, Aap and Xlp^/Xn for e — 0.05 in the case of EG without a strong nonlinear term. 
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0.005 


-0.598881 


1.66316 


-2.60733 


0.01 


0.340322 


-0.15760 


-0.98090 


0.02 


0.259368 


-0.22086 


-1.62996 


0.05 


0.119439 


-0.11385 


-1.86873 


0.1 


0.061951 


-0.06598 


-1.93656 


0.2 


0.031501 


-0.03117 


-1.96879 



Table 5: The values of A„, Aap and A^^/A„ for e = 0.1 in the case of EG without a strong nonhnear term. 



6 Discussion 

In Ref. the authors argue that there is a pressure difference at plate 2, namely the pressure in one side 
of the plate is different to that on the other side. They call this new pressure which acts on the central plate 
the "flux induced osmosis" (FIO). We consider the existence of FIO proposed by ^B] identifying AP as the 
pressure difference defined in Sec. |S| 

In Ref. [El the following criteria are stated: 

1. AP > regardless of the sign of Qy. 

2. Pyy is a function of the noncquilibrium quantities: Ti, the nonequilibrium steady heat flow Qy, and 
it is related to the equilibrium quantity Pgq as long as the noncquilibrium and equilibrium temperature at 
both sides of plate 2 coincide, 

"(0) _ (dPyy\ ^23^ 



where Ti is the thermodynamic temperature of plate 2. 

In Ref. 14^, for a system of hard spheres and of maxwellian particles — which obey Boltzmann's equation 
and which obey the BGK equation jl8| — using the Chapman-Enskog method, it is shown that criterion 1 in 
Ref. J5^ is valid but criterion 2 is not valid. 

On the other hand, for the hard disk's system, from our present scheme based on Enskog's equation we 
obtain that 

1. Criterion 1 in Ref. is not obeyed in the case of NLBG: AP is always negative independent of the 
sign of Qy. 

2. In the case of LEG, ((221) is vahd. However, substitution of ([TTIl and into leads to X^^ /Xn = -2. 
This is correct only in the case of NLBG. Hence, criterion 2 in Ref. 15 is not satisfied either. 

In the formulation of SST it is assumed that the number of particles and the size of the system is infinite 
but that the number density of the system is finite 15 . This condition implies that the terms 0{S^) in 
the collision terms can be neglected. Table |21 still indicates that under such condition the osmotic pressure 
difference is not always positive. Especially when the system is extremely close to equilibrium, Eq. (|21|) 
implies that AP is always negative. 

However when the system is extremely close to equilibrium, Eq. (|22|l implies that AP is always positive. 
This result coincides with the hard sphere case in Ref. |14| . 

Furthermore, the condition under which there is no heat flux at the porous wall is 

/°° /-oo j-oo pO 

-OO Jo ^ J —oo J -oo ^ 



(24) 



Using Eq. Q yields 



9y(0) = J^[neq-n(0)]. (25) 
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The above condition is not satisfied without introducing a difference between the temperature of the gas 
and the porous wall except in the LBG case. Therefore, except in the LBG case, it must be difficult to 
maintain the equilibrium state between plates 1 and 2 even if the heat conductivity of plate 2 is extremely 
high. Here, as a rough simplification, let us introduce the temperature of the gas in contact with plate 2, 
Tg^Ti. Substituting 5 = inY^qa.^, ^ anA using the no-mass-fiux 

condition given by an equation similar to yields, 



AP = 7i(0) 



1 - 



2JT„ 







+ 









Ql 



256n(0)' 



(26) 



As Ti < Tg < T2 it follows that 1 < Tg < 1 + e, and putting qy{yY = = in the right side of Eq. il'^) . not 
only for NLBG hut also for LBG, AP < 0. The behavior of the pressure difference changes qualitatively even 
if we restrict the analysis to Boltzmann's regime. This implies that the estimation of AP is a very delicate 
problem. Even if one can prepare the walls which satisfy Eq. H23|) and estimate the pressure difference in 
Enskog's regime, it is difficult to know what physical meaning lies behind such case. 



7 Conclusions 

We have analyzed a simple noncquilibrium steady state system inspired by T5^ . Our study refers only to 
a hard disk system and analyze in great detail its behavior using our extended hydrodynamic equations (2| 
using various approximations. Since we obtain that the osmotic pressure difference is negative in many cases 
for which Eq. (|23|l is not satisfied we cannot agree with 

We have assumed that the pores in plate 2 are small enough and we have not considered the problem 
about refiections on the wall at all. As we point out in the last part of Sec. the boundary (contact) 
condition is very delicate. We recognize that a more sophisticated analysis is necessary. 

However, in cases when strong nonlincarities can be neglected and the system is quite close to equilibrium — 
so that higher order terms in e do not contribute — then the osmotic pressure is positive. This implies the 
possibility of the existence of FIO. 

In addition, in the full-paper 15 , the authors point out that Eq. (|23() is directly related to the condition 
at the wall and this condition is essential to construct the formalism of SST in a complete form which gets a 
new noncquilibrium extensive quantity which determines the degree of noncquilibrium. Therefore, to clarify 
the problem, the measure of the pressure difference is done only for the case of a wall obeying Eq. H23|l . 
Within this context, they still recognize the results of ^1] as implying the existence FIO. Hence we also 
think that it is worth estimating the pressure difference starting from Eq. (|23|) . If one were to analyze 
the problem in such a way then the boundary (contact) condition would have to be reconsidered to solve 
the kinetic equation. In other words, one would have to evaluate AP under the rather complex conditions 
required by kinetic theory that would lead to satisfy Eq. H23f) . This has not been done. 

The SST formalism is quite interesting and the present study has only put to test the possible existence 
of FIO. 

Finally, we briefiy comment about EIT: Extended irreversible thermodynamics |19[ 1201 1^ . For an ideal 
gas, Refs. [201 1^ studied a problem quite similar to the one in the present article. In Ref. ^50" the authors 
estimated the pressure difference without considering a special wall. They assumed that the direction of the 
heat flow is parallel to the interface and their results are very interesting. Furthermore, in the authors 
estimated the difference between the pressures which are parallel and perpendicular to the heat flow. 

For a hard disk, appling LBG to the systems of |j2(ll21j . it is easily possible to get similar equations to 
® and to see that the pressure difference predicted by [201 is positive and that the difference predicted by 
PT] is zero. In both cases, of course, it is totally unnecessary to use conditions Q and Q. 
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